function ret=yz_rotation_sqrt(alpha,theta,beta)
 mat=yz_rotation_form(alpha,theta,beta);
 if(theta==0)
  ret=[alpha/2 theta beta/2];
  return;
 endif
 ret=my_sqrt(mat);
endfunction

